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ABSTRACT 

Context. Accretion disks in AGN should be subject to the same type of instability as in cataclysmic variables (CVs) or in low-mass 
X-ray binaries (LMXBs), which leads to dwarf nova and soft X-ray transient outbursts. It has been suggested that this thermal/viscous 
instability can account for the long term variability of AGNs. 

Aims. We test this assertion by presenting a systematic study of the application of the disk instability model (DIM) to AGNs. 
Methods. We are using the adaptative grid numerical code we have developed in the context of CVs, enabling us to fully resolve the 
radial structure of the disk. 

Results. We show that, because in AGN disks the Mach numbers are very large, the heating and cooling fronts are so narrow that 
they cannot be resolved by the numerical codes that have been used until now. In addition, these fronts propagate on time scales much 
shorter than the viscous time. As a result, a sequence of heating and cooling fronts propagate back and forth in the disk, leading only 
to small variations of the accretion rate onto the black hole, with short quiescent states occurring for very low mass transfer rates only. 
Truncation of the inner part of the disk by e.g. an ADAF does not alter this result, but enables longer quiescent states. Finally we 
discuss the effects of irradiation by the central X-ray source, and show that, even for extremely high irradiation efficiencies, outbursts 
are not a natural outcome of the model. 
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1. Introduction 

Accretion disks are found in a large variety of astronomical ob- 
jects, from young stars to active galactic nuclei. Among these, 
close binaries have deserved special attention, because they are 
nearby, and vary on short timescales that enable time-dependent 
studies of their light curve. In particular, a number of these sys- 
tems show large outbursts, as for example dwarf novae, which 
are a subclass of cataclysmic variables in which a low-mass com- 
panion transfers mass onto a white dwarf; these systems undergo 
outbursts lasting at least a few days during w hich the ir bright- 
ness increases by several magnitudes (see e.g. lWarrie r 1995lfor 
a review). The outbursts are belie ved to be due to a thermal- 
visco us accretion disk instability dMever & Mever-Hofmeisteil 
Il98ll) which arises when the disk effective temperature becomes 
of order of <10 4 K, enough for hydrogen to become partiall y ion- 
ized a nd opacities to depend strongly on temperature (see lLasotal 
2001 for a review of the model). Similarly, soft X-ray transients, 
which are a subclass of low-mass X-ray binaries in which the 
compact object is either a black hole or a neutron star also show 
outbursts, their amplitude being larger and the time scales longer 
than for dwarf novae. The ionization instability of the accretion 
disk is also thought to be cause of the outbursts, the difference 
with dwarf novae being due to the difference in the mass of the 
compact object (and thus in the depth of the gravitational poten- 
tial well) and to the effect of illumination of the disk , much more 
impor tant in the case of X-ray binaries (see e.g. iDubus et al.l 
120011) . 



It was realized long ago (Lin & Shields 1986) that the same 
instability could be present in accretion disks around AGNs; it 
was found that, at radii ~ 10 15 ~ 16 cm where the effective tem- 
perature is indeed of a few thousand degrees, the disk should be 
unstable. For the parameters of AGNs, the implied timescales 
are of order of 10 4 - 10 7 yr, making impossible the direct obser- 
vation of the instability, but predicting that in many systems the 
disk should not be in viscous equilibri um and that many AGNs 
should be in a quiescent stat e (see ISiemiginowska et al.l li 996; 
Siemiginowsk a & Elvislll 99% . It was also immediately realized 
that, as in dwarf-novae, the character of putative AGN outbursts 
strongl y depends on assumptions o ne makes about the disk vis- 
cosity dMineshige & Shi elds 1990). However, while in the case 
of dwarf-novae one is guided by the observed outburst proper- 
ties when fixing the viscosity prescription, in the case of AGN 
it is not even clear that outbursts are present as the variability of 
these objects could be due just to mass-supply variations. This 
state of affair gave rise to various, more or less arbitrary, pre- 
scriptions of ho w viscosity varies (or not) with the state of the 
accre t ion flow dMineshige & Shields! 1 199(it iMenou & Ouataertl 
120011; Uaniuk et alJ l2004)". In addition, results of numerical cal- 
culations of AGN outbursts were marr ed by the insufficien t res- 
olution of grids used. As showed bv lHameurv et alj dl998l) low 
grid resolution often leads to unreliable results^ 
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1 Maver & Pringle (2006) mischievously remark in this context that 
"mathematical convergence does not necessarily imply more accurate 
modeling of physical reality". While this might be true it is clear that 
the lack of convergence of a mathematical model makes it useless for 
physical applications. 
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The aim of the present article is the systematic analysis of the 
application of the DIM in the context of AGN disks. Instabilities 
other than the thermal-viscous instability may exist in AGN 
disks (beyond the MRJ instability which i s thought to be the 
source of viscosity iBalbus & Hawlevl 1 199 lb and in particular 
the gravitational instability that arises when self-gravity exceeds 
the c o mbined action o f pressure and Coriolis forces (Toomre 
Il964t ISafronovlll960l) : conditions for the onset of this insta- 
bility are met a t large distances from the black hole (see e.g. 
Shlo smanll990 ). The outcome of this instability in the AGN case 
is most probably the fragmentation of the acc retion disk (see e.g. 
Gam miel200 fl iGoodmanteOOl lRafikovl2"007h since in the AGN 
case, t he cooling time is likely to be short. iDuschl & Britschl 
(2006) suggested that the gravitational instability might instead 
be a source of turbulence, which could be the case if the non- 
linear development of the instability does not lead to fragmen- 
tation, not a likely outcome in the AGN case as mentioned 
above. Other local or global i nstabilities may arise, such as the 
Lightman-Eardley instability (Lig htman & Eardlevll 19741) . but it 
is far beyond the scope of this paper to discuss them all, and we 
consider parameters such as these instabilities do not occur. 



2. Vertical disk structure 

We recall here the vertical-structure equations adapted to AGN 
parameters. We consider here only the case where the viscosity v 
is proportio nal to the gas pressure (not th e total pressure, in order 
to avoid the Lightman & Eardley ( 1974) instability). The vertical 
structure of an a disk in which the viscosity v is assumed to be 
proportional to th e gas pressure is given by the standard disk 
equations (see e.g. Frank et al. 2002 and references therein): 



dP 

dz 



~P8z 



— = 2p 
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= V, 



dlnT 
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where P — P g + f ra d, p and T are the total (gas plus radi- 
ation) pressure, density and temperature respectively, g is the 
surface column density between vertical coordinates -z and +z, 
g z = Q^z the vertical component of gravity, Qk being the 
Keplerian angular frequency, F z the vertical energy flux and V 
the temperature gradient of the structure. This is generally radia- 
tive, with V = V ra d, given by: 



V,, 



rad 



kPF z 



(5) 



When the radiative gradient is superadiabatic, V is convective 
(V = V conv ). The convective gradient is calcula ted in the mix- 
ing le ngth approximation, in the same way as in lHameurv et alj 
(1998), with a mixing length taken as H m \ = a m \Hp, where Hp 
is the pressure scale height: 



Hp 



pg z + (Pp)Vm K ' 



(6) 



which ensures that Hp is smaller than the vertical scale height of 
the disk. Here, we use a m i = 1.5. 



Note that we have neglected the disk self gravity. This ap- 
proximation is valid as long as the ratio of self gravity to that of 
the central object is small: 



< 1 



(7) 



If this not the case, the disk is gravitationally unstable, which, as 
mentioned in the introduction, is likely to lead to fragmentation 
if the cooling time is short enough, or may significantly change 
the angul ar momentum transport by introducing non local t erms 
(see e.g. iLin & Pringlel Il987t iBaibus & Papaloizoul [l999l) . In 
both cases, the thermal-viscous instability can no longer apply 
(in the first case for obvious reasons, and in the second one 
because non-local effects cannot be approximated by viscosity, 
which is local); in our calculations we always make sure that the 
condition (|7]) is fulfilled. 

The parameter a s g is an effective viscosity, equal to the stan- 
dard viscosity coefficient a when the disk is in thermal equi- 
librium, but which also accounts for the time-dependent terms 
which are assumed to be also proportional to the pressure (see 
iHameursT et al. 1998 for a detailed discussion). 

The equation of state of matter is interpolated from the tables 
of iFontaine et al.l (1 19771) : in the low temperature regime (below 
2000 K), which is not covered by these tables, Sa ha equ ations 
are solved iteratively, as described by iPaczvnskil d 19691). The 
Rosseland mean opacities are taken from | Cox & Tabor! (Il976l) 
above 10,000 K, and from lAlexanderl (1 1975b below (more mod- 
ern opacities intro duce changes that are not important in the 
present context, see lLasota et al.ll2008l) . 

The boundary conditions are q — and F z = at the 
disk midplane, and q — S at the surface. The standard pho- 
tospheric condition KP g = 2/3g z has to be slightly modified, 
as (1) radiation pressure can de dominant, and (2) g z can vary 
in the photosphere. Integrating the vertical hydrostatic equilib- 
rium equation, and using the Eddington approximation leading 
to T 4 (t) = 3/47^(2/3 + t) where t is the optical depth, one 
obtains 



2 1 1 

3 \ KpZ 



(8) 



The term 1/Kpz is of the order of the relative thickness of the 
photosphere relative to the total disk thickness; it is usually of 
little importance except when the disk luminosity is close to its 
local Eddington limit, in which case the photosphere can be quite 
extended. 

The thermal equilibrium corresponds to Q + - Q~, where 
Q + and Q are the surface heating and cooling rates respectively 
(see Eq. ( TBl i below). Figure[T]presents two examples of thermal 
equilibrium curves in the E - r e jf plane, showing the charac- 
teristic S shape. Also plotted are the conditions g s /g c = 1 and 
h/r = 0.1. As can be seen, self gravity becomes important at 
radii l arger than 1 - 2 X 10 16 c m, i n agreement with the find- 
ings of lCannizzo & Reifj d 19921) and lCannizzol d!992l) . The con- 
dition that self-gravity be small can be quite sev ere; we note 
for ex ample that in several of the simulations by iJaniuk et al.l 
(l2004i) this condition is not fulfilled and the corresponding re- 
sults are therefore invalidated. The thin disk approximation con- 
dition (h/r <sc 1) is usually less stringent; it may however break 
for high accretion rates, in which case radiation pressure gradi- 
ent almost balances vertical gravity in a significant fraction of 
the disk vertical extent. 
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Fig. 1. Examples of S curves in the 2 - T eft - plane, for M = 1O 8 M and r = 10 15 cm (left) and for r = 2 x 10 16 cm (right). In both 
cases, S curves obtained for a — 0. 1 and 0.01 are shown. The dotted curve correspond to ratio of self to central gravity equals 1, and 
the dashed curve to h/r - 0.1. Only regions below the dashed curved and above the dotted one are allowed. For r - 10 15 cm, these 
limits lay outside the portion of the 2 - T e g plane shown here. 



The values 2 m ; n and 2 max that are the minimimun (resp. max- 
imum) values of 2 on the upper (resp. lower) branches of the S 
curve can be fitted by: 

,1.04 
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cm 
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cm 



-2 



(9) 



(10) 



and by: 

I.„=3. S 510V»»(^I-)"V,.» g 

where M 8 = M/(1O 8 M ), M being the black hole mass; 
the corresponding effective temperatures are r ef f(2 m ; n ) = 
4300(r/10 15 cm)-° 12 K and r eff (2 max ) = 3300(r/10 15 cm)-° 12 
K respectively. Note that these are independent from a, as ex- 
pected, and that their radial dependence is quite weak. As com- 
pared to disks around stellar mass objects, the surface densities 
are much higher and hence the effective temperature somewhat 
smaller at the turning points of the S -curve (the upper stable so- 
lution ends at 3000 - 4000 K instead of 7000 - 8000 K), even 
though the corresponding mid-plane temperatures are quite sim- 
ilar. 

It should also be noted that, because we are restricted to a re- 
gion where self gravitation is small, the disk extension, as mea- 
sured by the ratio ri n /r out is not very large: for the case of a 10 8 
M black hole, this is of order of 100, i.e. comparable to disks in 
CVs, but much smaller than for LMXBs. It is also worth noting 
that the disk thickness H : 



H 

r 



Vk' 



(11) 



where c s and \\ are the sound and Kepler velocities, is small in 
regions where the thermal-viscous instability can propagate. As 
compared to the CV or LMXB case, c s is unchanged because the 
central temperature at the turning points of the S curve are simi- 
lar, and of the order of 10 6 cm s _1 , but is much larger since in 
the AGN case we are restricted to regions close to the black hole 
as mentioned above; here, we have H/r ~ 10~ 3 or smaller. As 
the width of the heating and cooling fronts are of order of a few 
times H dMenou et al.l fl999), this can be a source of numerical 
problems. In particular, these fronts have been by far unresolved 
in all previous studies, casting some doubt on their results. 
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Fig. 2. Vertical structure of the accretion disk, for M = 10 8 M o , 
r = 10 I5 cm, 2 = 4.3 x 10 4 g cnr 2 , a = 0.01 and T c = 1400 K 
(thick curve) and 1570 K (thin curve). The line is dashed when 
energy transport is radiative. These two cases correspond to the 
upper and lower stable part of the small wiggle shown if Fig.Q] 



2. 1 . Critical points and the viscosity prescription 

It should be noted that for some choices of parameters the equi- 
librium curves show two "wiggles" (see e.g. the case a = 0.1 at 
r = 10 15 cm, FigJTJ on the lower branch. This also happens in 
accre tion disks around s tellar mass black holes, but, in contrast 
with iJaniuk et al.l ((2004), we do not find that this is always the 
case for AGNs. This discrepancy can be due to a difference in 
the treatment of convection, or to different opacities. The small 
wiggle at low temperature is not related to a strong change in 
the opacities, but instead to a strong change in the adiabatic gra- 
dient when molecular hydrogen becomes partially dissociated, 
as is shown by Fig. [2] Two vertical structures, corresponding to 
the same r, 2, and a, but two different effective temperatures on 
the upper and lower stable branch of the small wiggle of Fig. 
[1] differ essentially by a much stronger temperature gradient in 
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the disk midplane: in one case, V a( j ~ 0.10, while in the other 
V a d ~ 0.40. This effect occurs only if molecular hydrogen be- 
comes partially dissociated in the convective zone; since, for the 
corresponding temperatures and densities the opacities are rela- 
tively low, there are cases where the transition between molec- 
ular and atomic hydrogen occurs in a radiative zone, in which 
case no wiggle is found. 

In the standard dwarf-nova model, it is assumed that the a- 
parameter changes rapidly when the disk temperature reaches 
the ionization instability; this is required for the amplitude of the 
modeled outburst to be comparable to the observed one. It is of- 
ten stated that the physical reason for such a change is the change 
in the ionization parameter of the gas, hence a is assumed to re- 
main constant when transiting this secondary wiggle. This also 
seems to be a reasonable hypothesis in the AGN case, and does 
not require the physics of accretion disk to be different in dif- 
ferent environments, even though the te mperatures and den sities 
are similar. In the following contrary to lJaniuk et al.l (2004), we 
assume therefore that the critical E max of the cold stable branch 
corresponds to the ionization instability and that the lower wig- 
gle is not associated with a change in a. T his point is of im- 
portance, since as shown by iHameurvl (|2002|) . the shape of the 
resulting S curve and hence the outcome of the model is by far 
dominated by the change in the viscosity parameter a. 

Finally, one should remark that whatever arguments are 
used to justify the change in a, i.e. the use of an a co \<i and 
an ofhot ~ (4 — 10)ar co id the real reason is the necessity to 
produce the required out burst amplitude. It has been argued 
dGammie & Menoulll998l) that the difference between viscosi- 
ties in the high and low (quiescent) states of dwarf-nova disks 
is due to the "decay" of the MRI mechanism that is supposed to 
be the source of turbulence in accretion disks (Balbus & Hawleyl 
1 1 99 lb . In the environment of AGN disks, the MRI is supposed 
to be operating also in c old disks, which was use d to argue that 
in this case o-hot ~ a^oid dMenou & Ouataert 2 0011) . 

However, as noted by Steven Balbus (private communica- 
tion), because of the fact that numerical simulations treat the tur- 
bulent dynamics of disks at a level far beyond anything that can 
be approached with strictly analytic techniques, there has been a 
tendency to grant simulations a level of certainty that they do not 
yet merit. A careful treatment of realistic energetics still remains 
beyond the capabilities of current codes, and even simple poly- 
tropic shearing box calculations need to be run at much higher 
resolutions and for much longer times than were once thought 
necessary. 

Therefore the values of critical Re ynolds numbers d e duced 
from numerical simulat ions only dGammie & Menoul 1 19981; 
Menou & Ouataert 2001) are highly uncertain and we opted for 
using the standard dwarf-nova DIM also in AGNs. 



where v r the radial velocity in the disk, j = (GM\r) 1 ^ 2 is the 
specific angular momentum of material at radius r in the disk, 
Q.K = {GM\ /r 3 ) 1 / 2 is the Keplerian angular velocity 

T he energy conservati on equation is taken as (see lCannizzol 

119931; IHameurv et ail 19981 for details): 



dT c 
~dt 



2(q + -q- + j) p c id(rv r ) ar, 



p c Cp r dr 1 d 



V 1 . (14) 



where P c and p c are the midplane pressure and density, and Q + 
and Q are the surface heating and cooling rates respectively. 
They are usually taken as Q + = (9/8)vEQ 2 c and Q - o~T^, 
r e fj being the effective temperature. The term J accounts for the 
radial energy flux carried by viscous processes: 



/ = l/rd/dr(rF e ). 



(15) 



where F e is the flux carried in eddies with characteristic velocity 
v e and size / e , is: 



dT c 
CpY,v e ——l e 
Or 



3 „ r dT c 

-vCpE— 
2 Or 



(16) 



These are identical to the equations of a disk in a binary system, 
except that there are no tidal torques and no tidal dissipation. 

The inner boundary condition is also unchanged from the 
binary case: 



vE = at r — r„ 



(17) 



where r m is the radius of the inner edge of the disk, and can be 
larger than the radius of the innermost stable orbit if the disk 
is truncated by the formation of an ADAF, in wh ich case n n is a 
given function of the mass accretion rate (see e.g. Hameury et al. 
1997). As we use ln(2) as a variable, the 2 = boundary condi- 
tion is not applicable. Instead, we take: 



vS = 1.1S„ 



at r — n r 



(18) 



so that this allows both for the disk to be in the hot or cold state 
(see below for a more detailed discussion on the effect of using 
this boundary condition). 

The outer boundary condition is more problematic, as the 
disk extends to large distances where all the usual approxima- 
tions are invalid (thin disk, neglect of self gravity, ...). We instead 
assume that at some distance r oul ~ 10 16 cm for M = 10 8 M Q , 
the mass transfer rate is given and constant. This approximation 
is valid provided that the heating front does not reach this outer 
radius. 

The heat equation (TT4l > requires two additio nal boundary 
conditions; as discussed in IHameurv et al.l d!998l) . these are of 
little importance, and we take J — at r — r m and r = r out . 



3. Disk evolution 

3.1. Basic equations 

The standard equations for mass and angular momentum con- 
servation in a geometrically thin accretion disk can be written 

as: 



as 13 

^7 = T-(^'r) 

ot r or 
and 



i a / 3 



ax i a - . , - 



(12) 



(13) 



3.2. Results 

Figure[3]shows an example of the evolution of the accretion disk. 
We have considered here a 10 8 M black hole accreting at 10 24 
g s _1 , about one hundredth of the Eddington limit: 



^Edd 4nGMm v 



Edd 



C 2 T] 



0- th CT] 



1.4 x 10 26 M 8 g s" 1 



(19) 



where 77 ~ 0.1 is the efficiency of accretion, an d crxh is the 
Thom pson cross section. The mass is the same as in J aniuk et al.l 
d2004l) . but the mass transfer rate is lower by a factor ~ 10 in 
order to avoid (or try to avoid) the instability to occur in the self 
gravitating part of the accretion disk. We have assumed a vary- 
ing a between ah = 0.2 and a c = 0.04. We have also taken the 
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Fig. 3. Time evolution of an accretion disk with the following 
parameters: black hole mass: 10 8 M , inner and outer radius: 
10 14 and 10 16 cm respectively, and mean mass transfer rate: 10 24 
g s . Top panel: visual magnitude, intermediate panel: accre- 
tion rate onto the black hole, lower panel: radius at which the 
transition between the hot and cold regimes takes place. 
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Fig. 4. Time evolution of an accretion disk with the following 
parameters: black hole mass: 10 8 M , inner and outer radius: 
10 14 and 10 16 cm respectively, and mean mass transfer rate: 2 
10 22 g s _1 . Top panel: visual magnitude, lower panel: accretion 
rate onto the black hole. 
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Fig. 5. Details of the outburst shown in Fig. [4] Top panel: vi- 
sual magnitude, intermediate panel: accretion rate onto the black 
hole, lower panel: radius at which the transition between the hot 
and cold regimes takes place. The red-dashed line is the semi- 
analytic value of the minimum transition radius given by Eq.l23l 



outer disk radius to be 10 16 cm, as the disk becomes self gravi- 
tating at larger distances. As can be seen, the disk can never be 
brought completely to the cold regime; as a consequence, rel- 
atively low amplitude oscillations are seen in the visual mag- 
nitude and in the accretion rate onto the black hole. This sit- 
uation is reminiscent of what happens in the case of soft X- 
ray transients when no disk truncation or irradiation is assumed 
(Men ou et all l2000t jDubus et al.ll2001h . or in the case of sym- 
biotic stars dDuschllll986l) : a cooling front is reflected at some 
radius much larger than the disk inner radius and as a result 
a heating front starts propagating outwards but it cannot quite 
reach the outer disk edge and a cooling front forms again. Such 
reflections occur when the surface density S behind the cooling 
front reaches Z max which triggers a new instability. The result- 
ing heating front propagates outwards until the post-front den- 
sity reaches £ m j n . Then a new cooling front starts going down 
the disk. There is however a significant difference in that the 
short time oscillations do not result in oscillations of the mass 
accretion rate on the same time scale. M fluctuates only on the 
longer time scale of the front oscillation pattern. The basic rea- 
son for this is that the front propagates at approximately a times 
the sound speed, i.e. on a time scale 



r r 

thtmt = = T^th; 

ac s h 



(20) 



where f t h is the thermal time scale. ff ront is shorter than the vis- 
cous time f v ; sc = (r/h) 2 t± by a factor r/h, i.e. by several or- 
ders of magnitude. The cooling front therefore propagates so 
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Fig. 6. Radial structure of the disk. The blue solid line represents 
the surface density, the red dashed one the central temperature. 
The dotted lines are the critical 2 m ; n and E max . See text for de- 
tails. 



rapidly that the surface density at smaller radii does not change; 
to a first approximation, it cannot propagate in regions where 
£ > 2 max (a c ). In the CV case, ff ront is shorter than f v j sc , but not 
by such a large amount and strong gradients in the disk make the 
effective viscous time comparable to the front propagation time. 

It must also be noted that the front occasionally reaches the 
outer disk edge; then the outer boundary condition which im- 
poses in particular that there is no outward mass flow is not valid, 
so that the correct sequence is probably different; the back and 
forth propagation of heating fronts on a short time scales is how- 
ever a firm prediction of the model. 

For lower mass transfer rates, the outer part of the disk can 
remain on the cold, stable branch, in which case the front prop- 
agation is restricted to the innermost parts of the disk. Figures 
@] and [5] show the evolution of disk with the same parameters as 
in Fig. [3] but with a mass transfer rate of 2 10 22 g s _1 , 5 10~ 3 
times the Eddington limit. As can be seen in Fig. [5] heating and 
cooling fronts propagate in a restricted fraction of the accretion 
disk. They do not reach radii larger than that at which the disk 
can sit on the stable cool branch, given the externally imposed 
mass transfer rate. They also do not reach the innermost regions 
where the surface density always remain high enough for the 
disk to be stable on the hot branch, except when entering a qui- 
escence periods, which happens when the disk finally empties 
on a much longer viscous time. The active phase lasts for 3 10 5 
yr in the case presented here, with more than 400 consecutive 
oscillations. Note that these are not random, but show relatively 
regular sequences of decreasing oscillations that are interrupted 
by an oscillation with a larger amplitude, clearly visible in Fig. 
[5] Note also some sort of hierarchical structure of these oscilla- 
tions. Figure [6] shows the radial structu re of the disk durin g the 
oscillating phase (compare wih Fig. 2 of Du bus et alj 2001). The 
semi-stable inner and outer regions are clearly visible; the un- 
stable zone is in fact divided in two regions: an inner unstable 
one, and an outer marginally stable one, where 2 « 2 m i n , result- 
ing from the successive passage of heating fronts that die at radii 
decreasing with time; a leftover of the death of these fronts is 
the little wiggle in S that gets smoothed with time as a result of 
diffusion, or when a heat front is able to reach this region. Note 
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Fig. 7. Instability of the inner disk edge when the boundary con- 
dition 2 = is used. The top panel shows the position of the 
transition radius as a function of time, the lower panel the mass 
accretion rate onto the black hole. The parameters are those of 
Fig. |4] These large amplitude are unphysical, since mass flows 
into the black hole from a region much larger than the width of 
the zone in which these fluctuations occur. 



also the spike in the unstable region, that carries a small amount 
of mass that will cause the small wiggles in the marginally stable 
region. 

3.3. Minimum radius reached by cooling fronts 

The minimum radius reached by the cooling front can be deter- 
mined by noting that the front propagates down to a point where 
2 = £max(«cX an d that the innermost parts of the disk are in 
quasi viscous equilibrium, which means that the surface density 
is determined by the accretion rate, almost constant in this hot 
inner region. This is equivalent to stating that, at the reflection 
point, the dissipation rate Q + is: 



3GMM 4 
= g7rr 3 / = crr e 4 ff (S max ,hot) 



(21) 



where / — 1 - (r/r; n ) -1 ^ 2 . Note that T e g is calculated on the hot 
branch, and is not given by the analytic fits obtained in Section 2. 
An examination of Fig. Q]shows that r e ff(X max , hot) is about 3.2 
times that of the turning point on the cool branch, r e ff(S max , cold) 
for a c /ah — 0.1. As Q + is proportional to vS, and hence propor- 
tional to a, one can guess that: 



r e 4 ff (£ max ,hot) = 3.2 4 -^r e 4 ff (i max ,coid) 



(22) 



which is also a very good approximation even for a c /ah = 1, as 
can be seen from Fig.Q] Now, from the fits of r 4 ff (E max , hot) we 
have: 



r=1.7x 10 1 



M a c 
10 Zi g s 1 ah 



(23) 



Figure|6]shows that there are two points where 2 crosses the 
Smax line, and Eq.|23]has indeed two solutions, one for which / 
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Fig. 8. Long term evolution on an AGN accretion disk when 
ri n can vary as a result of e.g. evaporation or the formation of 
an ADAF. The system alternates between active phase in which 
heating and cooling fronts propagate back and forth in the disk 
and quiescent phases lasting about one million years. The top 
panel: visual magnitude, intermediate panel: accretion rate onto 
the black hole; lower panel: inner disk radius. 



is small, and another one in which f - 1. The first one corre- 
sponds to the transition between the very inner disk, where 2 is 
vanishingly small because of the boundary condition and there- 
fore the cool branch solution applies, and nearby regions where 
£ is large enough for the hot solution to apply. This will be dis- 
cussed in the next section. The second one corresponds to the 
radius at which the cooling front is reflected into a heating front. 

From Eq. [23] it appears that for low enough M or large 
enough r; n , the cooling front can reach the inner radius, in which 
case the system will enter a quiescence phase. More precisely, 
this happens when Eq.[23]has no solution. Simple algebra shows 
that the critical M is: 



M = 5xl0 20 (-^f 5 ^ 
\10 14 cm/ a c 



g s- 1 



and the corresponding critical radius is r — 1.44/?j, 
/ = 1/6. It is interesting to rescale this relation as: 



M 

M EdA 



2 . 7 10-^Mf (r m /3r s ) 2 - 5 



(24) 
at which 

(25) 



which shows that, if the disk is not truncated, low states will be 
found only for low mass transfer rates, whatever the black hole 
mass. 

The critical rate given by Eq.[25]refers to the accretion rate 
onto the black hole, and not to the mass transfer rate. For high 
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Fig. 9. The outburst shown in Fig. [8] The two upper panels are 
the same as in Fig. [8] the lower panels shows the transition radius 
between the inner hot disk and the cooler outer parts 



mass transfers rates, both are almost equal, as seen above (see 
e.g. Fig. [3); for lower values of M tlans f , they may differ by up to 
one order of magnitude, as shown for example in Fig. [4] where 
the accretion rate at maximum is 10 times larger than the mass 
supply rate to the disk. In this case, M ultimately falls below 
the critical value given by Eq. [25] and the disk enters a quies- 
cent state. The duration of this state is short, however, of order 
of the duration of the outburst state, as the average mass accre- 
tion rate during the active state is ~ 5.7 10 22 g s _1 , i.e. not very 
different from the steady mass transfer rate (note that the duty cy- 
cle expected for an outburst with an average accretion rate 2.85 
times larger than the transfer rate is 0.35, very close to the value 
given by the simulation (0.30), showing that the disk is almost 
relaxed). We therefore expect that outbursts exist only for low 
mass transfer rates, that these outbursts are weak - never reach- 
ing anything close to the Eddington limit - and that the duty 
cycle cannot be large. 

3.4. Innermost disk instability 

The very inner parts of the disk, where the density is very low 
because of the inner boundary condition should therefore be on 
the cold branch; the transition between this cold region and more 
distant, hotter regions should also be unstable. Indeed, when one 
assumes that the inner boundary condition is not S = l.lE m ; n 
at r — r; n , but is smaller than £ m j n , oscillations are found. For 
the sake of completeness, we show in Fig. [7] the effect of such 
oscillations in such a case. Cooling/heating fronts propagate in a 
very restricted region, whose radial extent is larger than the ver- 
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tical scale height so that the thin disk approximation is still valid, 
but much presumably smaller than the zone from which matter 
flows into the black hole. Note also that this region is so small 
that the total disk luminosity remains constant. These oscilla- 
tions are possible only when the width of heating/cooling fronts 
is smaller than the width of the region over which S catches the 
boundary condition S = 0, which is of the order of one to a few 
percents of r m (see Fig. |6), otherwise fronts would simply not 
exist. This is quite possible in the AGN case, because the fronts 
are so narrow and contrasts with the CV or LMXB case where 
the reverse is true and the condition E = does not have such 
an effect. There oscillations are most probably not physical, be- 
cause one assumes that (i) there is absolutely no torque at the 
inner disk edge, and (ii) that matter is lost from the disk only 
at r — n n ; it is very likely that the mechanisms leading to ac- 
cretion at the inner disk edge (e.g. evaporation, ...) will smooth 
oscillations there. In order to avoid these, and in order to ease 
the numerical computations, we have assumed that S is not van- 
ishingly small at the inner disk edge, but that instead it is very 
slightly larger than £ m j n . 

3.5. Disk truncation 

Disk truncation could be a solution to the absence of large out- 
bursts; this was found to ba a n essential ingredient of the soft 
X-ray transient model (see e.g. iMenou et alj|200& iDubus et alj 
120011) . Truncation can be the result of the formation of an 
advection-domin ated accretion flow (ADAF), or of one of its 
variants (see e.g. lNaravan & McClintockll2008t iKato et al.ll 19981 
for reviews of the ADAF); the important feature being that the 
flows becomes hot, geometrically thick, and optically thin close 
to the black hole. In order for the outburst cycle to be modified, 
one needs the disk not to extend down to the innermost stable 
orbit, but instead be truncated at a radius comparable to the min- 
imum radius reached by the cooling front. The inner disk radius 
will then depend on the mass accretion rate onto the black hole. 
Many prescription can be derived; what really matters is whether 
Eq. d23l can be satisfied or not, the details of the variations of r; n 
as a function of M being of no importance. 

Figures [8] and [9] show an example in which r m - 2 10 14 
(M/10 23 g s -1 ) 1 / 2 . As can be seen, quiescent states are found, 
as well as active state which are not very bright though - only 
brighter than the active states for non truncated disks by a factor 
~ 2, resulting in duty cycles which also differ by factors ~ 2. 



3.6. Disk irradiation 

Disk irra diation plays an es s ential role in soft X -ray transients 
(see e.g. Ivan Paradiis|[l996t IDubus etalJl200ll and references 
therein), and sometimes in CVs, see e.g. "jHameurv et al.lll999b ; 
it could also play an important role in the AGN context. We fol- 
low here the same procedu re as in the case of irradiated disks 
in SXTs (Dubus et al. 2001). We assume that the irradiation flux 
Fin- onto the disk is given by: 



Fit = &IZ = C 



Mc 1 
4jtR 2 



(26) 



where C is a co nstant; in the case o f SXTs, C = 5 x 10~ 4 has 
been adopted bv lDubus et alj d200ll) (note that we include here 
in C the efficiency of mass to energy conversion). The ratio of 
the irradiating flux to the viscous flux is then given by 
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Fig. 10. Time evolution of an irradiated disk, with parameters 
identical to those of Fig. [4J here C = 0.01. Top panel: V mag- 
nitude; intermediate panel: accretion rate onto the black hole; 
lower panel: transition radius between the hot, inner disk and 
the cool outer disk. 
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Fig. 11. Example of the disk radial structure in the irradiated 
case. The blue solid line is the surface density, the red dashed 
line is the central temperature. The blue dotted lines represent 
Smin and £ max , and the effects of irradiation are clearly visible ar 
radii < 10 15 cm; at radii smaller than 3 10 14 cm, the S shape of the 
cooling curve vanishes and the disk is stable 



4 c M acc r 
3 M r s 



(27) 
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which clearly shows, because we are interested in regions much 
closer to the black hole than in the case of SXTs, that C must be 
large if irradiation is to have any effect at all: in order to affect the 
central temperature of the disk, Fi n /F v i sc must exceed the optical 
thickness of the disk (in the radiative case). This is in principle 
possible because the X-ray emitting region could have a complex 
geometry, as e.g. a corona above a cooler disk, in which case the 
irradiation flux can be large; it is however very unlikely that it 
could exceed 0.01, since the fraction of X-rays absorbed below 
the photosphere is at most 10%. In the following, we consider 
the case C = 0.01, which corresponds to a maximally irradiated 
disk. One should also note that in this case, the X-ray luminos- 
ity would be linked in a complex way to the local properties of 
the accretion flow, and would not be proportional to the accre- 
tion rate onto the black hole M acc only. Note also that Eq. (|27T i 
assumes steady state, and neglects the r m /r terms in the energy 
dissipation equation, which can be significant close to the disk 
inner edge. 

We have then calculated a grid of vertical disk structures in 
order to determine the effective temperature as a function of 2, 
T c , T m , as described in Section 2, with a modified boundary con- 
dition at the disk surface 

F z = cr(T 4 - 7£) (28) 

As in the case of SXTs, the effect of irradiation is a stabilization 
of the disk when the irradiation temperature is large enough, typ- 
ically larger than ~ 10 4 K. 

A reasonably good fit to the effective temperature at E max is 
given by: 

- ™U=) <"(lo^r) 

in a situation where the disk is in viscous equilibrium and hence 
Jin- is directly given by M, valid for irradiation temperatures 
larger than about 4000K, but less than 10 4 K for which the disk 
becomes stable. 

Figure[10]shows the time evolution of an irradiated disk with 
the same parameters as in Fig. [4] apart from the irradiation factor 
C set to 0.01. As can be seen, even in this maximally irradiated 
disk, the time evolution is not very different from that of unirra- 
diated disk. There is still a succession of rapid oscillations of the 
luminosity, with a heating/cooling front propagating back and 
forth; the main difference is here that the disk enters into a qui- 
escent phase more rapidly than in the unirradiated case, but one 
should note that the initial structure was not exactly the same in 
both cases, and that because of the huge computing time required 
to follow the disk oscillations, a relaxed state can not be attained 
in practice. However, the radial disk structure obtained at the end 
of active phases in the irradiated and unirradiated case, with an 
outer disk on the cold stable branch, and most of the disk having 
2 = S m ; n ; this similarity, and the fact that the unirradiated disk 
was almost relaxed makes us confident that here also the disk is 
close to relaxation. 

Figure QT| shows the radial structure of the disk at time t = 
3.5 10 4 yr, when 3/4 of the first outburst have elapsed. It clearly 
shows the impact of irradiation on the innermost part of the disk, 
which is due both to a large irradiation temperature and to the 
decreasing viscous dissipation close to the inner disk edge (the 
/ factor). 

The conclusion that quiescent states are possible only for low 
mass transfer rates is very general; an analysis similar to that 



described above in the non irradiated case leads to the conclusion 
that the disk can enter into quiescence only if the mass accretion 
rate is less than 

*»- = 4 *^(l0^Mgf V <30) 

accurate to within a factor 2 when compared to the results of the 
numerical simulations. Although the dependence on Mg, r m , and 
a are different from the unirradiated case, the numerical value 
of this critical rate is not changed to the point where one could 
obtain large amplitude outbursts, during which the Eddington 
limit would be attained or approached. 

4. Conclusion 

We have shown that the accretion disks in active galactic nuclei 
can indeed be subject to the same thermal-viscous instability as 
in dwarf novae and soft X-ray transients, but the outcome of this 
instability is very different. This contrasts with previous findings 
that large amplitude outbursts reaching the Eddington limit were 
possible, and the reason for this discrepancy is the insufficient 
spatial resolution of the numerical codes that have been used 
to model the disk. In AGNs, the disk opening angle is much less 
than in DNs or SXTs, because the Keplerian velocity is not small 
compared to the speed of light, whereas the sound speed is, by 
construction, the same in both cases. This results in very thin 
transition fronts which are quite difficult to follow numerically. 

We do however predict time variations of the AGN lumi- 
nosity by a few magnitudes on time scales ranging from a few 
thousand (the propagation time of a thermal front in the disk) 
to a few million years (the typical quiescent/outburst time), and 
these oscillations are enhanced by a possible truncation of the 
innermost parts of the disk. However, because of their too small 
amplitudes and duty cycles, these variations cannot explain the 
statistical properties of quasar and AGN luminosity distribution. 

One should also stress that, for high mass transfer rates, 
the transition front can reach regions in the disk where self- 
gravitation becomes important, and there the assumption of a 
homogeneous disk becomes quite questionable, as a result of the 
development of a gravitational instability that is likely to result 
in the fragmentation of the disk. 
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